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Abstract 

We investigate the quantum fluctuations of the fields produced in sum-frequency (SF) 
generation from light initially in the photon number squeezed state. It is found that, to the 
fourth power term, the output SF light is sub-Poissonian whereas the quantum fluctuations 
of the input beams increase. Quantum anticorrelation also exists in SF generation. 


1 Introduction 

In recent years squeezed states of light have been successfully produced in several laboratories 
and certain applications also demonstrated . 1 In future applications it may be desirable to change 
the frequency of a light beam already in the squeezed state, for example by frequency doubling, 
parametric down conversion or four wave mixing. Three wave processes are preferable since the 
second order nonlinear susceptibility is much larger than the third order susceptibility. Second 
harmonic generation from a pump beam initially in a squeezed state has been discussed 2,3 as 
well as SF/DF generation and degenerate parametric down conversion of quadrature squeezed 
light . 3,4 In this paper we shall consider SF generation from input fields initially in the photon 
number squeezed or sub-Poissonian state. 

2 Sum-frequency generation 

As shown in Fig. 1 we take both input fields, of frequency and Cl> 2 > and the sum-frequency 
(J 3 field — UJ 2 + Wi) to be in a collinear geometry. The input beams are incident on a 
nonlinear medium such as a crystal with an effective second order nonlinear coefficient x** The 
field Hamiltonian for ideal SF generation ignoring losses is 

H = hujia\ai + ^u; 2 a 2 a 2 + huj^a^a^ + &X f ( a i a 2<4 + h.c.) (l) 

where a u a 2 and a 3 are the annihilation operators for the two input beams and the resultant 
SF beam, respectively. 
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FIG. 1 Sum-frequency generation in a nonlinear crystal. 

We introduce the slowly varying operators 

e,(t) = a.j (2) 

where the subscript j = 1, 2 or 3 refers to the two input or SF waves, respectively, as throughout 
our paper. For propagation along the z-axis at a velocity u in a lossless medium we make the 
conversion z — ut. Since u> 3 = uq + wj and substituting x = X* /** we deduce the equations of 
motion for the slowly varying operators from Eq. (1) to be 

* 4 “ = 

4 ¥ d = -'ViW'jW. ( 3 ) 

iJd = -xcWcW. 

In the short path approximation for x z 1 we assume the solution of Eq. (3) to take the form 
Cy (z) = c ; (0) + {x z )Bj + {xtfCj + (x z )* D i + iX z Y E } + ‘ * * ( 4 ) 

where c,(0) denotes the incident fundamental field at t = 0. It may be seen that c,(z) still 
satisfies the commutation relation 

[«(*).«}(*)] = */• ( 5 ) 

Substituting into Eq. (3) we obtain the solutions for c^z ), c 2 (z) and c 3 (z). 

To discuss the photon number fluctuations of the quantized fields we consider the variance 
(An 2 (z)) or the Fano factor 

< 6) 

<W)> 

where n,(z) = c}(z)c,(z), and <An?(z)) = {nj(z)> - (n^z)) 2 . To obtain the above expressions 
we need to solve for the factors c\c k , c* 2 c* and c* 3 cj, where c k = c*(0) and k = 1 or 2, then find 
their expectation values. If the input fields are in a photon number squeezed state we adopt the 
formalism given by Kitagawa and Yamamoto 5 , that is, 

Ci» = e i ' 1kr ' k a k + fc, e* = -**?*«*«•■♦* . (7) 


Here n k = a\a k is the number operator and ^ k , T) k are constants chosen to maximise the squeezing 
for a given input field intensity proportional to |a*| 2 . Using Eqs. (7) and after some tedious 
calculations we may obtain the expectation values (c^c t ) and higher terms with respect to the 
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eigenstates of a*, noting that Cs(0) is a vacuum field. The Fano factors may thus be computed 
numerically for explicit values of a*, 7* and 77*. It is evident that the expressions for the Fano 
factors of the two input beams are identical if the subscripts k = 1 and 2 are interchanged. 

Another quantity that reflects the quantum fluctuations of the light fields is the quantum 
correlation between the sum of the input photon numbers, 711(2) + 712(2), and the SF photon 
number 713(2). That is, if we define (71(2)) = (711(2) + 712(2) + 713(2)), then 

(Ati 2 (2)) = <A(n a + n 2 ) 2 (2)> + (A Tig (2) ) + 2 [<( 7 Z X + 71 2 )( 2 )ti s ( 2 )) - <( 7 i! + n 2 ) (2)) (713(2))] , ( 8 ) 


and we define the Fano factor for the total sum of the photon numbers as 


0(2) = 


(An*(2)) 

<*( 2 )> • 


( 9 ) 


The Fano factors for different pairs of the two input beam intensities are plotted against the 
effective interaction length \z in Figs. 2 to 4 . It can be seen that a s becomes less than 1 in 
all three cases, that is, the SF produced becomes sub-Poissonian. Moeover, as the interaction 
length increases the photon number fluctuations are reduced even further but at the expense of 
the input fields. 




FIG. 2 Fano factors of the input and SF beams versus effective interaction length 
for sub-Poissonian inputs of equal intensity. Subscripts 1,2 and 3 denote the two 
inputs and the SF beam, respectively, a is the Fano factor for the total photon 
number of the three beams. 

As shown in Fig. 2 , when the two input beams are of the same intensity, that is, |<*i| J = ] 02 1 2 
and assuming that their initial statistical properties are identical, then and a 2 a re always 
the same, as may be expected. However, when one frequency e.g. has a greater number of 
photons than the other, its Fano factor is less affected (see Fig. 3 ) since there is a greater number 
of unconverted photons left. The Fano factor a for the sum of all three intensities also incre ases 
with interaction length, but it is always better i.e. less than that for the individual beams. This 
can be understood from the correlation between the three beams, since the creation of one SF 
photon is always associated with the annihilation of two input photons, one from each input 
beam. 
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FIG. 3 Fano factors of the input and SF beams versus effective interaction length 
for sub-Poissonian inputs of different intensities. Subscripts 1,2 and 3 denote the 
two inputs and the SF beam, respectively, a is the Fano factor for the total photon 
number of the three beams. 

In Fig. 4 we have taken one input beam u> 2 to be sub-Poissonian whilst the other Wj is in a 
coherent state, with both having the same intensity. It can be seen that the curves of a 2 and a 3 
behave similarly to those in the previous figure, but ai remains constant at 1. This implies that 
the coherent input beam remains coherent throughout the interaction, with only the originally 
squeezed beam suffering increased quantum fluctuations. The physical reason for this is yet to 
be fully understood. 



FIG. 4 Fano factors of the input and SF beams versus effective interaction length 
for inputs of equal intensities but different statistical nature. Subscripts 1,2 and 3 
denote the two inputs and the DF beam, with cuj in a coherent state and u 2 in a 
sub-Poissonian state, a is the Fano factor for the total photon number of the three 
beams. 
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3 Conclusion 


We have investigated the photon number fluctuations in SF generation with sub-Poissonian light 
as input and small effective interaction length. We find that the generated SF beam becomes sub- 
Poissonian and the quantum fluctuations decrease as the effective interaction length increases, 
but at the expense of the two input beams which become noisier. The Fano factor of the sum 
of the intensities of the three beams increases also with interaction length, but it will always 
be less than that of the other beams. This can be understood on the grounds of the quantum 
correlation between the input and SF fields. When one of the input beams is coherent and the 
other sub-Poissonian, the coherent beam remains coherent throughout the interaction, whilst 
the other becomes noisier. The generated SF beam, however, still becomes sub-Poissonian. 
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